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ABSTRACT 


A  finite  elsment  analysis  of  the  end  notch  flexure  specimen  for  mode  II 
interlaminar  fracture  toughness  measurement  was  conducted.  The  effect  of 
friction  between  the  crack  faces  and  large  deflection  on  the  evaluation  of 
Gjjq  from  this  specimen  were  investigated.  This  analysis  showed  that  the 
friction  effect  is  negligible  when  crack  length  is  greater  than  one-fourth 
of  length  of  the  specimen.  Further*  geometric  nonlinear  analysis  is  required 
to  evaluate  Gjj^  when  deflection  is  greater  than  the  thickness  of  specimen. 


INTRODUCTION 


The  delamination  of  laminated  fiber  reinforced  composites  Is  of 
current  interest,  particularly  with  regard  to  their  durability  and 
damage  tolerance.  It  is  the  most  commonly  observed  failure  mode 
In  composites .  Fracture  mechanics  has  been  found  to  be  a  useful 
tool  for  understanding  delamination  [I],  Most  of  the  work  reported 
In  this  area  has  been  concerned  with  the  delamination  in  the  opening 
mode,  I.e.  mode  I.  Delamination  under  shear  sliding  mode  I.e.  mode  II 
and  combined  mode  l-ll  situations  are  also  the  possible  fracture  modes 
of  potentially  greater  significance  in  composite  structures.  Several 
studies  on  mixed-mode  interlaminar  fracture  are  also  reported  [1-41. 

Recently  attention  has  been  given  to  mode  11  shear  delamination 
[5-9].  Several  investigators  have  used  end  notch  flexure  (ENF) 
specimen  to  measure  shear  mode  strain  energy  release  rate,  C,,p,  for 
delamination  growth  [6-9].  This  specimen  is  a  beam  with  a  crack 
(delamination)  located  on  the  neutral  plane  at  one  end  and  is 
subjected  to  the  three-point  bending.  This  specimen  was  initially 
utilized  by  Barrett  and  Foschi  to  measure  Mode  11  critical-stress- 
intensity  factor,  for  a  strongly  anisotropic  fibrous  material, 

wood  [10J.  Currently  an  ASTM  task  group  is  conducting  a  round 
robin  test  program  using  an  end  notch  flexure  specimen  as  c  prelude 
to  developing  an  ASTM  standard  for  measurement  of  mode  II  interlaminar 
fracture  toughness  of  composites.  In  this  study,  the  end  notch  flexure 
specimen  was  analyzed  using  the  finite  element  method  for  understanding 
Its  behavior  and  performance. 


FINITE  ELEMENT  ANALYSIS 

The  end  notch  flexure  specimen  geometry  is  shown  in  Fig.  t. 

Two  geometries  of  the  specimen  were  analyzed;  t1)  total  length, 

2L  =  7.62  cm  and  (2)  2L  =  10.16  cm.  For  both  geometries,  the 
thickness  (2h)  of  the  specimen  was  taken  equal  to  0.30  cm  which 
is  nominal  thickness  of  a  24  ply  composite  panel.  These  dimensions 
of  L  and  h  were  selected  since  they  were  employed  in  experimental 
works  16,9].  The  width  of  specimen  was  2.54  cm.  Only  one  material 
system  was  analyzed;  unidirectional  graphite/epoxy.  The  material 
properties  used  in  finite  element  analysis  are  presented  In  Table  1. 

A  two-dimensional  finite  element  program  called  GAMNAS  [11] 
was  used.  A  typical  finite  element  mesh,  shown  in  Fig.  2  consisted 
of  about  1000  isoparametric  four-node  elements  and  had  about  2400 
degrees  of  freedom.  The  element  size  In  the  vicinity  of  the  crack 
was  0.02  x  0.02  mm  In  order  to  evaluate  accurately  the  strain  energy 
release  rate  Cjj.  This  size  was  selected  by  previous  experience  112} 
as  well  as  by  conducting  the  convergence  test  on  C||  calculation  with 
mesh  refinement.  The  error  in  C||  obtained  In  the  present  study  is 
estimated  to  be  less  than  ±2%.  The  strain  energy  release  rate,  G||, 

In  the  analysis  was  computed  using  a  virtual  crack  closure  technique 
[13J.  Plane-strain  condition  was  assumed  in  the  finite  element  analysis. 

The  analysis  indicated  that  mode  II  deformation  Is  achieved  at 
the  crack  tip  and  that  the  accompanying  mode  I  deformation  causes 
closure  and  overlapping  of  the  opposite  faces  of  the  crack.  This  is 
not  possible  physically.  To  prevent  the  overlapping  of  crack  faces, 
the  nodal  coupling  technique,  available  with  GAMNAS  program  111], 
was  used.  For  this  purpose,  the  multipoint  constraint  was  applied 
at  corresponding  nodes  to  have  the  same  displac  jments  normal  to  the 
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crack  faces.  This  resulted  in  pure  mode  II  condition,  i.e.  C|  =  0. 
When  the  crack  faces  are  pressed  together  as  In  this  specimen, 
friction  effects  will  occur.  The  friction  effects  on  Cjj  calculation 
will  be  rather  difficult  to  assess  in  the  practical  situation,  however 
a  reasonable  estimate  of  its  role  can  be  investigated  with  the  finite 
element  analysis.  The  friction  effect  was  modeled  in  the  following 
manner.  The  analysis  was  conducted  intiaiiy  with  the  frictionless 
crack  faces  where  displacements  normal  to  the  crack  faces  were 
matched  at  corresponding  nodes.  The  normal' reactions  at  the  crack 
faces  were  obtained  from  this  first  step  of  calculation.  The  friction 
forces  at  the  nodes  on  the  crack  faces  were  calculated  for  the  assumed 
coefficient  of  friction.  The  finite  element  analysis  of  the  specimen 
was  conducted  with  these  computed  friction  forces  applied  at  the  crack 
faces.  This  second  step  of  analysis  provided  the  new  normal  forces 
at  the  crack  faces  and  corresponding  friction  forces.  The  analysis 
was  repeated  till  the  applied  friction  forces  matched  with  the  computed 
normal  reactions.  This  required  about  two  or  three  iterations. 
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RESULTS  AND  DISCUSSIONS 


Comparison  Between  Theoretical  and  FEM  C,.  Values 


The  measured  compliance,  c,  versus  crack  length,  a,  relationship 
of  a  test  specimen,  in  general,  provides  the  most  accurate  evaluation 
of  fracture  toughness.  However,  it  requires  several  measurements.  On 
the  other  hand,  fracture  toughness  can  be  measured  from  a  single 
measurement  of  critical  load  if  the  necessary  expression  for  G  calculation 
is  available.  The  experimental  dc/da  values  were  found  to  be  subject  to 
large  errors  due  to  small  variation  in  compliance  with  crack  length  for 
small  crack  length  values  in  the  ENF  specimen  [71.  Therefore,  the 
following  expression  for  mode  II  strain  energy  release  rate,  Cjj  for 
this  specimen  was  developed  by  Russell  [7] : 

9P^a^ 

^11  ~  2  3 

where  P,  a,  b,  h  are  defined  in  Fig.  1  and  E^^  is  the  Young's  modulus 
in  longitudinal  direction.  The  C||,  as  expressed  above,  is  based  on  a 
linear  beam  tiieory  where  shear  deformation  is  neglected  and  the  curvature 

1  _  d^y/dx^ 

^  [1  +  (dy/dx)^j^^^ 

is  approximated  with  1 /Rti  d^y/dx^.  The  C||  values  obtained  from  FEM 
analysis  and  from  Eq.(l)  for  both  geometries  of  the  specimens;  (1) 

2L  =  7.62  cm  and  (2)  2L  =  10.16  cm,  are  compared  for  the  normalized 
crack  length,  a/L,  ranging  from  0.125  to  0.9.  The  normalized  crack 
length  commonly  used  for  measurement  of  with  this  specimen  was 
0.5  [6,7,9}.  Figure  3  shows  that  C||^  from  a  linear  beam  theory. 


I.e.  from  Eq.(1),  is  within  8%  and  6%  of  its  counterpart  from  FEM 
analysis  for  a/L  =  0.5  for  specimens  with  2L  =  7.62  cm  and  2L  =  10.16  cm, 
respectively.  Further.  Fig.  3  shows  that  deviation  between  Cjj  from 
Eq.(1)  and  FEM  analysis  decreases  with  longer  crack  lengths  and 
increases  with  shorter  crack  lengths.  This  comparison  shows  that  the 
Cue  ®^*^***  ^  evaluated  from  the  end  notch  flexure  specimen  for 
a/L  i  0.5  If  Eq.(1)  is  used. 

Effect  of  Overhang  Length 

The  effect  of  overhang  length  beyond  the  support  was  also 
investigated  with  finite  element  analysis.  Specimens  with  overhang 
lengths  of  6.35,  12.7  and  25.4  mm  on  each  side  were  analyzed.  The 
C||  values  obtained  for  these  cases  were  within  ±1%  of  each  other. 

Thus,  the  overhang  portion  of  this  specimen  has  no  effect  in  the 
evaluation  of  C||^. 

Friction  Effect 

As  mentioned  previously,  the  friction  was  also  modeled  in  FEM 
analysis.  The  most  difficult  task  is  to  evaluate  the  actual  coefficient 
of  friction,  u,  between  two  pressed  crack  faces  of  ENF  specimen. 

This  was  beyond  the  scope  of  the  present  study.  However,  a  simple 
experiment  was  performed  to  estimate  the  coefficient  of  friction  between 
two  unidirectional  graphite/epoxy  laminates.  This  experiment  involved 
the  measurement  of  the  force  to  Initiate  sliding  between  the  two 
laminates  with  a  known  compressive  load.  The  coefficient  of  friction,  ji, 
was  found  to  be  about  0.35.  Therefore,  the  finite  element  analysis 
was  conducted  with  the  three  values  of  p  =  0.1,  0.3  and  0.5  to  under¬ 
stand  the  role  of  friction  in  a  qualitative  manner.  Results  of  this 
analysis  are  shown  in  Fig.  4.  Several  observations  can  be  made  from 
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this  analysis.  For  crack  length  commonly  used  in  experimental 
investigation,  i.e.  a/L  =  0.5,  the  effect  of  friction  is  negligible 
for  all  practical  purposes  if  u  is  about  0.3.  And,  this  effect  is 
less  in  the  ENF  specimen  with  2L  =  10.16  cm  in  comparison  to 
specimen  with  2L  =  7.62  cm.  The  friction  effect  shows  a  pronounced 
effect  for  shorter  crack  lengths  and  is  negligible  for  longer  crack 
lengths.  This  leads  to  an  apparent  conclusion  that  longer  crack 
length  should  be  employed  for  evaluation,  but  this  may  cause 

the  large  deflection  of  specimen  during  testing  and  henceforth  it 
may  require  the  geometric  nonlinear  analysis  as  discussed  later  on. 

The  friction  forces  between  the  crack  faces  were  primarily 
in  the  region  above  the  support  pin.  Since  it  is  always  desirable 
to  have  the  friction  force  as  small  as  possible.  It  might, therefore, 
be  preferable  to  use  the  shim,  or  rollers  between  the  crack  face 
just  above  the  support  pin  to  reduce  the  coefficient  of  friction. 

Barrett  and  Foschi  [10]  suggested  the  use  of  roller  between  crack 
faces  for  measurement  of  K||^  of  wood  from  this  type  of  specimen. 

However,  it  may  not  be  convenient  in  composite  specimens  due  to 
their  small  thicknesses. 

Geometric  Nonlinear  Analysis 

To  investigate  the  effect  of  finite  rotation  due  to  large 
deflection,  the  geometric  nonlinear  analysis  of  the  ENF  specimen 
was  conducted  with  the  CAMNAS  finite  element  program  [llj.  This 
analysis  was  conducted  without  friction.  Figure  5  shows  the  comparison 
of  C||  obtained  from  linear  and  geometric  nonlinear  analysis  for  both  geo¬ 
metries  of  specimens  for  crack  length  a  =  0.5L.  Linear  and  nonlinear  solution 


are  in  agreement  with  each  other  up  to  deflection,  A ,  equal  to  the 
thIcKness  of  specimen,  2h.  For  larger  deflection,  I.e.  A  i  2h,  the 
Cji  would  be  overestimated  from  a  linear  analysis  or  Eq.(l).  This 
overestimation  depends  on  the  deflection  as  shown  In  Fig.  5. 
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CONCLUDiNG  REMARKS 

A  finite  element  analysis  of  the  end  notch  flexure  specimen 
for  mode  II  interlaminar  fracture  toughness  measurement  was 
conducted.  The  effect  of  friction  between  the  crack  faces  and 
large  deflection  on  the  evaluation  of  from  this  specimen  were 

Investigated.  Two  geometries  of  the  specimen,  and  one  material 
system  was  studied.  This  study  showed  that  this  specimen  should 
be  preferably  used  for  Cj|^  measurement  with  crack  length, 
a  ^  0.5L  and  deflection,  A  <  2h. 
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iilSfSlMMmuS: 


TABLE  1  -  CRAPHITE/EPOXY  MATERIAL  PROPERTIES® 


Modulus**,  CPa 

Poisson's  Ratio** 

^11 

^2 

S2 

''12 

'*23 

131.0 

13.0 

6.4 

0.34 

0.3S 

*•  ^3  ®  ^2'  '’l3  “  ''l2"  ®13  "  ^2 


b:  The  subscripts  1,  2  and  3  correspond  to  the  longitudinal,  transverse, 
and  thickness  directions,  respectively,  of  a  unidirectional  ply. 


function  of  normalized  crack  length. 

•  4 


s  a  function  of  normalized  crack  length. 


as  a  function  of  normalized  crack  length. 


s  a  function  of  normalized  crdck  length. 
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